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Abstract 

The maximal value of thermal conductivity K max of the perfect wurzite GaN crystal con- 
taining isotopes of natural abundance is estimated. Our upper limit of k = 4800 W/Km 
at T max = 32 K is smaller than calculated by Liu and Balandin k = 6000 W/Km and 
higher than obtained by Slack etal k = 3750 W/Km. The phenomenological dependence 
k oc T~ L43 obtained by Mion et al for the temperature interval 300 — 450 K is extended to 
200 — 300 K. For temperatures higher than T max the best fitting of our experimental data 
to Callaway's formula is obtained for Grueneisen's constant equal to 7 = 1.35. 
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In heterostructures GaN/AlGaN in electric field of the order of lOOkV/cm, the den- 
sity of dissipated power reaches 10 6 W/cm 2 H1I2I3L This means that gallium nitride 
based devices suffer from heating effects that significantly limit the performance of 
high electron mobility transistors and laser diodes thereby reducing the device's 
life and reliability. This means that the accurate knowledge of the gallium nitride 
thermal conductivity is crucial to understanding thermal effects. For this reason, 
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the investigation of the temperature behavior of thermal conductivity coefficient is 
an actual and reasonable problem. 

Recently four groups measured thermal conductivity of GaN in the wide tem- 
perature ranges. Luo et al flU have measured the thermal conductivity of both 
LEO and HVPE single crystals of GaN over the temperature range 60 — 300 K. 
Slack et al H5I6H performed measurements in the temperature interval 11 — 300 K. 
We measured thermal conductivity of massive high pressure growth dislocation- 
free crystalline sample containing about 10 20 per cm -3 point defects in the tem- 
perature interval 4.5 — 300 K and performed a preliminary analysis of temperature 
dependence of it [7J. We also measured heat capacity of GaN for temperatures 
between 5K and 300K, and as a result we established the value of Debye's tem- 
perature of GaN [8]. These measurements also allow us to obtain the temperature 
dependence of the phonon mean free path L eff . Measurements of thermal conduc- 
tivity for temperatures higher than 300 K were performed by Mion et al ||9l. 

Slack et al compared their experimental results to the Callaway formula consider- 
ing separately the contributions of longitudinal and transverse phonons. Therefore, 
their fitting procedure relays on six fitting parameters, namely two Debye's tem- 
peratures 6> L , 9 T Grueneisen's constants 7l, 7t and constants 6 L , 6 T related to the 
Umklapp processes (L and T stands for longitudinal and transverse phonons, re- 
spectively). The single set of their free adjustable parameters contained also the 
relaxation rate for the normal processes rjy . Liu et Balandin [10J considered ther- 
mal conductivity GaN. They argued that in the case of GaN specimens the normal 
processes do not influence the thermal conductivity because of strong scattering on 
point defects. For this reason we use Calaway's formula with k 2 = 0. The results 
of measurements of thermal conductivity on wurzite GaN samples H4I5I6I7H were 
analyzed by Kamatagi et al [11] employing Holland's and modified Callaway's 
models. A good representation of the temperature dependence of k was obtained 
with a single set of phonon-phonon scattering prameters. The purpose of this study 
is different than Kamatagi et al [fTTTl . Having in mind the needs of high-power elec- 
tronic devices manufactured from GaN we shall establish the upper limit of thermal 
conductivity for T > T max . We base our approach on the standard Callaway for- 
mula and use only two free adjustable parameters. 

In dielectric and semiconducting crystals heat is carried mainly by acoustic phonons. 
As a rule, the group velocities of such phonons are much greater than those of opti- 
cal phonons. Therefore, heat flow in semiconducting crystals is limited by various 
mechanisms of scattering of acoustic phonons. There are two principal kinds of 
phonon scattering in crystals, namely resistive and normal processes lfT2l . They 
are characterized by the corresponding relaxation times r R and r N . Normal pro- 
cesses do not change the total momentum of phonons. Such processes influence 
the heat flow only indirectly. The resulting phonons of higher energy are scattered 
more intensively on point defects and have greater probability of participating in 
Umklapp processes. In the perfect crystals (dislocation and chemical impurity-free 
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samples), the relevant intrinsic resistive scattering is due to three-phonon Umklapp 
with relaxation time t\j, the point defect scattering with characteristic time r P , and 
the boundary scattering r B . The characteristic time r P is determined mostly by the 
crystal natural isotope composition. 



(1) 



The total thermal conductivity k(T) can be written as 0~3 

K = Ki + K 2 , 

where K\(T) and k 2 {T) are defined by 
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where k p and H are Boltzmann's and Planck's constants respectively, v is the mean 
velocity of phonons, 6d is Debye's temperature, T is the ambient temperature, x 
frw/k B T, is the phonon frequency, tJ 1 
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When tn >> tr, the combined relaxation time r c m tr and the term k 2 © is 
negligibly small. The generally accepted expression for Umklapp phonon scattering 
rate is (cf. e.g. flU) 
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where 7 is Grueneisen's constant, M is the mass of GaN molecule, and is b is a 
fitting parameter. As a rule for crystalline solids, 2 < b < 3 lfl2l . 



For phonons scattered by isolated defects of mass different from that of the host in 
an otherwise perfect crystal, Tamura derived a scattering rate lfT4l 
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where V is the volume of the elementary cell, the parameter g is 
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with gr(cr) is the constant which express the strength of isotope effecs. It depends 
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on mass Mi(cr) and fraction fi(a) in the ith isotope of a atom as 
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fi{&)Mi(cr) is the average mass of a atom. 
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The boundary scattering is important for low temperatures and does not influence 
thermal conductivity in the intermediate and high temperatures. The related scat- 
tering rate 1 is 

(8) 
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where L cff is the effective phonon mean free path. For T <C 9 G it is of the order 
of the cross-sectional dimensions. A precise determination of L c q is crucial for 
correctly describing the k(T) curve in the low temperature range. 

As we have recently established, for GaN, the Debye temperature 9 = 365 K J8l, 
which is lower than used by other authors H5I61 101 1 511 6H . It is therefore obvious 
that one should revise the fitting procedure. We calculate g (Eq© accounting for 
natural composition of isotopes of Ga and N. Our calculations show, in agreement 
with ref. [0, that the contribution of nitrogen isotopes to g is ten times smaller than 
the contribution of isotopes of gallium. 

Since other defects always present in available specimens diminish the thermal 
conductivity coefficient, we obtain the upper limit of thermal conductivity of GaN 
specimens with the natural composition of isotopes. The Grueneisen parameter for 
hexagonal GaN was calculated by Lopuszyhski and Majewski lfi~7l . Their value 
7 = 1.28 is bigger that used by other authors H6I10I181I . For temperatures higher 
than T max (corresponding to /t max ), the best fitting is obtained for 7 = 1.35, which 
is close to the calculated value [fTTl . As a result of fitting the Callaway formula to 
our experimental results (with tn tr), we found b = 2.5. In TableQ]we collected 
the values of the above parameters. 
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9b [K] 


5.1 x 10 5 


1.35 


1.37 x 10" 4 


2.5 


365 



Table 1 

Single parameter set based on Callaway's model used to represent the temperature depen- 
dence of the thermal conductivity k 

Extension of temperature dependence calculated by us is in excellent agreement 
with k oc T -1 - 43 obtained in ref. (H in the temperature interval 300 — 450K for a 
specimen containing dislocations of the lowest density. This means that the phe- 
nomenological expression k = k(300K) x (T/300K) -1 ' 43 proposed in ref. |fT9l 
can be also used in temperature interval 200 — 450K. One should stress that this 
phenomenological expression differs from the results set forth in ref. [6l (cf. Fig. 
[I])- The observed numerical values of k and the temperature slope of k(T) indicates 
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that 3lu dynamic method of thermal conductivity measurements [9| and the axial 
stationary heat flow method measurements yield the same temperature depen- 
dence, hence these two experimental techniques should be considered as supple- 
mentary. 




Fig. 1. Thermal conductivity k(T) of GaN for T > T max . Open circles - our results Q, 
crosses are taken from ref. [9 ] for sample with the low dislocations density, filled squares - 
ref. 0. Solid line - the present calculations, dotted line - k oc T~ L22 ref. Q. 



In Fig. [2] we present the results of our calculations for 3 < T < 500 K and 
compare them with the results of calculations presented in refs. H6I9I10H . In the low 
temperature limit, where the boundary scattering is essential, the calculated k was 
fitted to our experimental results using the value of L cS = 0.08 cm. This value is 
close to the low temperature mean free path of acoustic phonons reported in ref. 

m. 

The mean value of velocity is calculated using the set of elastic constants reported 
by Polian et al [|20l . The obtained value of v is close to velocity v = 6 x 10 5 cm/s 
measured in our experiments on propagation of phonon pulses at helium tempera- 
tures 11211 

Our upper limit of k (K max = 4800 W/Km at T max = 32 K) is smaller than cal- 
culated by Liu and Balandin K max = 6000 W/Km IfTOl and higher than value 
K max = 3750 W/Km announced in ref. 0. 

In conclusion, for hexagonal GaN crystals containing isotope atoms and otherwise 
perfect crystal, we obtained an acceptable representation of all our data for temper- 
atures between 3 K and 500 K. Using Callaway's model with k 2 = (r N 3> r R ), 
we adjusted two free parameters (7 and b) in the combination of the two scattering 
mechanisms considered, with the scattering rates for isotopic and boundary scatter- 
ing given in Eqs. © - ®. The remaining parameters were calculated or deduced 
from independent experiments. We established that thermal conductivity for tem- 
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Fig. 2. Thermal conductivity of gallium nitride in the temperature interval 5 — 450K. Open 
circles - results for our best samples [7], thick solid line - present calculation, crosses - 
ref. [9], curve No. 1 - ref. [ 10], curve No. 2 - ref. (H. 

peratures greater than T max is determined by scattering by point mass defects and 
the Umklapp processes. In the temperature range 200 — 500 K we justified the 
empirical dependence k oc T~ L43 observed by Mion et al flU. 
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